Intermediate-Coupling Theory for the Spectral Weight of Spin Polaron 
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Within the intermediate-coupling theory, the quasiparticle weight Z of one hole injected in the 
undoped antiferromagnetic ground state is studied. We find that, for the hole located at the quasi- 
particle band minimun with momentum ko = (±5, ±5), Z is finite. By comparing the results 
obtained by the self-consistent Born approximation, we show that the intermediate-coupling theory 
for Z is appropriate only when J/t ~ 1.6. Finally, the reason why this approach fails in the small-J 
case will also be clarified. 



PACS Numbers: 75.10.Jm, 71.10.+X, 74.65. +n 

The problem for the motion of a hole in a 2D antifer- 
romagnet (AF) has received significant attention jjj es- 
pecially since the discovery of the copper oxide supercon- 
ductors, where superconductivity arises from the doping 
of holes in an antiferromagnetic insulator. The AF with 
one hole is also a highly non-trivial correlated electronic 
system, and is therefore of fundamental interest from a 
theoretical point of view. 

Intuitively one would expect that the presence of a hole 
in an AF leads to a distortion of the underlying spin con- 
figuration, in a similar way as a conduction electron in a 
polar crystal causes a deformation of the lattice. Indeed, 
under the assumptions that the AF is not completely 
destroied by doping a hole and that the low-energy exci- 
tations of the spin background are spin waves, one arrives 
at the following effective Hamiltonian H, which is 

reminiscent of Frohlich's polaron Hamiltonian B : 



H = H t + H Jt 

H.J = ^CJ q & q & q , 



(1) 



q 



Ht = 5f £( M q( k )/k- q /i4 + h - c -) 

V k,q 

Here 6 q and /k are the annihilation operators of the hole 
and the spin wave. a-> q = ujv^ is the spin wave exci- 
tation spectrum, where lu = Jz/2, f q = ^Jl — 7 q , and 

7 q = J^d e lq d /z with d the unit vectors to nearest neigh- 
bors and z being the coordination number ( z = 4 for 
a 2D square lattice ). a = 2t/J is the dimensionless 
coupling parameter ( therefore, the small-J limit means 
the strong-coupling limit of H ), N is the number of 
the lattice sites, and M q (k) = cosh(9 q 7k- q + sinh# q 7k 
is the coupling function between the hole and the spin 
wave, where cosh6> q = [(1 + v^) / {2v q )] 1 / 2 and sinh# q = 
— sgn(7 q )[(l — Vq) I (2t' q )] 1 / 2 are the Bogoliubov transfor- 
mation coefficients. Based on this similarity [pi, a hole 
in an AF may be viewed as a " spin polaron ", i.e., 
a hole dressed by a cloud of virtual spin-wave excita- 
tions of the antiferromagnetic spin background. From 



this observation, the intermediate-coupling treatment || 
of the Frohlich polaron problem is recently applied to 
the present spin-polaron problem by Barentzen Q. It 
is found that : (1) the intermediate-coupling results for 
the quasi-particle energy E{k) is in agreement with the 
dispersion curve obtained by means of a Green function 
Monte Carlo method ||; (2) the result for the band- 
width W is quite good for weak coupling ( J/t > 3 ), 
and is still reasonably good in the intermediate range ( 
0.4 ~ J/t < 3 ), where the deviation from the values ob- 
tained by the self-consistent Born approximation ( SCBA 
) H was about 10 - 20 %. Thus the intermediate-coupling 
theory may be appropriate for J/t ~ 0.4. 

One of the most controversial issues in the spin-polaron 
problem is whether a hole injected in the undoped ground 
state behaves like a quasiparticle [ jlO| ], or, equivalently, 
whether the quasiparticle weight ( or the wavefunction 
renormalization constant ) Z at the Fermi surface is finite 
under the dressing by the spin- wave excitations. There 
have already been many studies along this line (see Q for 
further references). However, most of the previous stud- 
ies have involved numerical calculations on small clusters 
( even the studies using the SCBA have to solve Dyson's 
equation numerically for small clusters ). Although nu- 
merical calculations on clusters show that the hole has a 
finite quasiparticle weight, there is still some uncertainty 
as to whether the quasiparticle weight vanishes or not in 
the thermodynamic limit ]il[| . 

In this report, we will study the quasiparticle weight 
within the intermediate-coupling theory, in which we can 
freely take the thermodynamic limit. Since, when a sin- 
gle hole is doped, the hole will locate at the quasiparticle 
band minimun with momentum ko = (±7£,±1|) we 
will confine ourselves to the spectral weight of the hole at 
momentum k . We find that the deviation from the re- 
sults of Refs. D is below 20 % only when J/t ~ 1.6. Thus 
the range of validity for Z is smaller than that for W. 
Moreover, we will point out that, although this approach 
is not plagued with the finite-size effect, due to its mean- 
field nature, the infrared behavior of the present system 
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may not be correctly described by this method even qual- 
itatively ! This may be the reason why the intermediate- 
coupling theory for the spin-polaron problem is not as 
successful as that for the lattice-polaron case. 

In order to compare with the results obtained by SCBA 
P|, rather than starting from the Hamiltonian used in 
Barentzen's paper ( Eq.(ll) of Ref. j7|, which is denoted 
as Hsar in the present report ), we take Eq.(|l]) as our 
starting point, where the Bogoliubov transformation has 
been taken, such that the unperturbated ground state is 
the vacuum state for the spinless fermion operators and 
the quantum Neel state, |0), with respect to the spin- 
wave operators 6 q . The relation between H and H Bar is 
simply 



H = V^H Bar V, 



(2) 



where V denotes the unitary operator of the Bogoliubov 
transformation. 

Following the procedure in Ref. 0, one first makes 
a change of coordinates to the rest frame of the mov- 
ing hole by the unitary operator ( i.e., the so-called Jost 
transformation ) |[12i 



c/ = E(^E^ ipRi )/^-p- 



1 N 



(3) 



where the translation operators for the bosons are defined 
by 



Ti = exp(-iR j; • ^2 q fo q fo q 



(4) 



Then, in order to further diagonalize all terms linear 
in the boson operators of the transformed Hamiltonian, 
U^HU, a displacement transformation W is employed: 



k 

^ k = cxp[-^^A q (k)(6 q -^)], 



(5) 
(6) 



where the unknown parameters A q (k) are determined by 
the mean- field equations via the variational principle, i.e., 



SE(k) 

SAq(k) 



= 0, 



(7) 



with the ground-state energy E(k) defined by the expec- 
tation value 



E(k) = (fl0\W^U^HUW\fl0). 



(8) 



It can be shown that, for a single hole, the Jost transfor- 
mation U and the Bogoliubov transformation V commute 
each other, then 



{fl0\W*U*HUW\fl,0) 
= (ft,0\W^U^H Bar VUW\fl0) 
= (fl,0\WWtfH Bar UVW\fl0), 



(9) 



which is just the expression of the expectation value ob- 
tained in Ref. JtJ . Thus the variational calculations in our 
case is completely the same as those in Ref. 0. There- 
fore, the mean-field equations lead to the following self- 
consistent equations: ( see Eqs.(57)-(58) of Ref. M ), 



M q 2 (k) 



a 2 o x - MiQg) 
lnF(k) = -^-F 2 (k)]T 



(^ q + 2afi(k)) 2 ' 



(10) 
(11) 



where F(k) and f2(k) are defined in terms of A q (k) as 

F(k)=exp(-l^A q (k)), (12 
q 

Q(k) = iF(k)^A q (k)M q (k). 



(13) 



Then the ground-state energy can be written in terms of 
F(k) and fi(k) : 



E(k) = -awfi(k)[l - 21nF(k)]. 



(14) 



Now we turn to the calculation of the spectral weight. 
The spectral weight of the hole at momentum k is 
defined as 



z k H<C'l<MI 2 , 



(15) 



where 



= |/k)0) is the ground-state eigenvector of 
the unperturbated Hamiltonian, Hj, within the one- hole 
subspace; while |^k) = UW\fl,0) is the corresponding 
variational ground-state eigenvector of the full Hamilto- 
nian, H. By Eqs.©-©, 



<°I^E T ^|o> 

i 

(0|W k |0) 

ex p(-^E A q ( k ))' 



(16) 



where 2^0) = |0) ( because 6 q |0> = ) and the Baker- 
Hausdorff formula is used in the last line of the derivation 
. From Eqs.@, ©, and @, one obtains 



Z k = exp(-i^A q (k)), 



= F(k). 



(17) 



Thus, by solving F(k) and f2(k) via the self-consistent 
equations, Eqs.(ll^) and (O), we get at the same time. 



2 



From now on, we confine ourselves to the case of 
k = ko. After numerically solving the self-consistent 
equations, we get Z = Z^ a as a function of J/t. The 
result is shown in Fig.l. If we compare our results ( solid 
line ) with those obtained by SCBA || ( open circles ), 
we realize that our spectral weight is fairly accurate for 
J/t > 2.5 ( i.e., the agreement is within 10 % ). However, 
the deviation from the results of Refs. |j] is below 20 % 

only when J/t ~ 1.6. Thus the range of validity for Z is 
smaller than that for W . 



FIG.1 
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FIG. 1. The spectral weight Z for one hole located at the 
quasiparticle band minimun with momentum ko = (±J,±J) 
as a function of J/t. The results of the intermediate-coupling 
theory are shown as the solid line, and open circles represent 
those obtained by SCBA §. 

In the following, we will clarify the reason why the 
intermediate-coupling theory for the spin-polaron prob- 
lem is not as successful as that for the lattice-polaron 
case. Notice that, after restricting to the subspace with 
one hole at momentum k, as shown in Eqs.(45c) and (47) 
of Ref . , one arrives at the effective Hamiltonian for the 
boson operators |l3| 

Hi;. ~ constant + e q (k)6 q & q , (18) 
q 

where e q (k) = ^ q + 2af2(k) can be considered as the 
renormalized energy spectrum of the spin waves in unit 
of lo in the rest frame of the hole at momentum k. 
Hence, for q — > 0, the renormalized energy spectrum 
e q (k) — * 2af2(k), which is nonvanishing for a given k; 
while the " bare " one, z^ q oc |q| — > 0. That is, a fi- 
nite gap, 2aJ7(k), is introduced into the excitation spec- 
trum of the spin waves by this variational approach, even 
though the original excitation spectrum is gapless due to 
the Goldstone theorem |L4| ! Thus the long-wavelength ( 
infrared ) properties of the system is altered qualitatively. 



This qualitative change in the excitation spectrum is not 
reasonable, because the correction to the spin wave en- 
ergies by a single hole in a macroscopic AF background 
should be proportional to 1/N, which would be negligible 
in the thermodynamic limit ! In the lattice-polaron case, 
although there is still a nonvanishing correction to the 
phonon energies §], there is no qualitative change in the 
phonon energies by this approach, because only the lon- 
gitudinal optical ( LO ) phonons are considered and the " 
bare " energy spectrum of these LO phonons can approx- 
imately be taken as a positive constant. As claimed by 
Anderson Jl0| , whether Z is zero or not greatly depends 
on the infrared behavior of the system. Since the infrared 
behavior may not be faithfully described by this varia- 
tional approach, it is reasonable that the intermediate- 
coupling theory may not predict the accurate value of Z 
in the spin-polaron case. 

For a further support of the above arguements, we 
present the results of the "induced spin gap", 2aCl(ko), 
as a function of J/t in Fig. 2. We find that, as J/t de- 
creases, the "induced spin gap" increases and invalidates 
the variational approach. 

FIG.2 
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FIG. 2. The " induced spin gap ", 2af2(ko), as a function 
of J/t. 

In conclusion, within the intermediate-coupling theory, 
we show that the quasiparticle weight Z is finite, and our 
results agree with those obtained by the self-consistent 
Born approximation |9| when J/t ~ 1.6. Because of the 
failure to describe correctly the infrared behavior of the 
system, this approach is not suitable for the study of the 
spectral weight of the spin-polaron case, especially when 
J/t is small. 
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